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We study the topological and fermionic vacuum structure of four-dimensional QCD on the lattice by means of 
correlators of fermionic observables and topological densities. We show the existence of strong local correlations 
between the topological charge density and the quark condensate, charge and chiral density. By analysis of 
individual gauge configurations, we visualize that instantons (antiinstantons) carry positive (negative) chirality, 
whereas the quark charge density fluctuates in sign within instantons. 



Over the last two decades several models have 
been developed to describe the basic properties 
of QCD, namely quark confinement and chiral 
symmetry breaking. The most popular are the 
dual superconductor, leading to confinement, and 
the instanton liquid model, which explains chiral 
symmetry breaking and solves the Ua(X) problem 
jlj . Both models rely on the existence of topolog- 
ical excitations, monopoles and instantons. In- 
stantons have integer topological charge Q which 
is related to the zero eigenvalues of the fermionic 
matrix with a gauge field configuration via the 
Atiyah-Singer index theorem [Q. Apart from this 
famous connection of topology and fermionic de- 
grees of freedom, here we attempt to systemati- 
cally shed light on the relationship between the 
sea-quark distribution and topology |§. We do 
this by studying correlators of topological densi- 
ties with fermionic observables of the form tpTtf) 
with T = 1,74,75. Those quantities are usually 
referred to as the quark condensate, quark charge 
density, and the chiral density. 

For the implementation of the topologi- 
cal charge on a Euclidian lattice we re- 
strict ourselves to the so-called field theo- 
retic definitions which approximate the topolog- 
ical charge density in the continuum, q(x) = 
JL-2^p° Tr (F ftv (x)F pa (xj\ . We used the pla- 
quette and the hypercube prescription. To get rid 
of quantum fluctuations and renormalization con- 
stants, we employed the Cabibbo-Marinari cool- 
ing method. Mathematically and numerically the 



local quark condensate xjjip(x) is a diagonal ele- 
ment of the inverse of the fermionic matrix of the 
QCD action. The other fermionic operators are 
obtained by inserting the Euclidian 74 and 75 ma- 
trices. We compute correlation functions between 
two observables 0\(x) and 02 (y) 



g(y-x) = (0 1 (x)0 2 (y)}-(0 1 )(0 2 ) 



(1) 
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and normalize them to the smallest lattice separa- 
tion d min , c(y-x) = g(y-x)/g(d min ). Since topo- 
logical objects with opposite sign are equally dis- 
tributed, we correlate the quark condensate with 
the square of the topological charge density, and 
similarly for the other quantities. 

Our simulations were performed for full SU(3) 
QCD on an 8 3 x 4 lattice with periodic bound- 
ary conditions. Applying a standard Metropolis 
algorithm has the advantage that tunneling be- 
tween sectors of different topological charges oc- 
curs at reasonable rates. Dynamical quarks in 
Kogut-Susskind discretization with rif = 3 fla- 
vors of degenerate mass m = 0.1 were taken 
into account using the pseudofermionic method. 
We performed runs in the confinement phase at 
(3 = 5.2. Measurements were taken on 2000 con- 
figurations separated by 50 sweeps. 

Figure |l| shows results for the correlation func- 
tions of Eq. (|j) with 0\ a local fermionic observ- 
able and O2 a topological density. All correla- 
tions exhibit an extension of several lattice spac- 
ings. Although cooling of quantum fields is nec- 
essary to extract topological structure, the cor- 
relations between the topological charge density 
and both the local chiral condensate and the ab- 
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Figure 1. Correlation functions of topological 
charge density with the quark condensate, quark 
charge and chiral density shown in the first, sec- 
ond and third row, respectively. 



solute value of the quark charge density are nearly 
cooling- independent. However, cooling (or some 
other kind of smoothing) is inevitable to obtain 
nontrivial correlations between the chiral density, 
Ci = ^Ibi'ix), and the topological charge den- 
sity. This can be expected since both quantities 
are correlated via the anomaly. 

We now turn to a direct visualization of 
fcrmionic densities and topological quantities on 
individual gauge fields rather than performing 
gauge averages. We persue this in the following 
to get insight into the local interplay of topol- 
ogy with the sea-quark distribution. By analyz- 
ing dozens of gluon and quark field configurations 
we obtained the following results. The topological 
charge is hidden in quantum fluctuations and be- 
comes visible by cooling of the gauge fields. For 
cooling steps no structure can be seen in q(x), the 
fcrmionic observables or the monopole currents, 
which does not mean the absence of correlations 
between them. After a few cooling steps clusters 
of nonzero topological charge density and quark 
fields are resolved. For more cooling steps both 
topological charge and quark fields begin to die 
out and eventually vanish. 

In Fig. [2] a typical topologically nontrivial con- 
figuration, consisting of an instanton and an an- 
tiinstanton, from SU(3) theory with dynamical 
quarks on the 8 3 x 4 lattice in the confinement 
phase is shown after 15 cooling steps for fixed time 
slices. We display the positive/negative topolog- 
ical charge density by white/black dots if the ab- 
solute value exceeds certain minimal fluctuations. 
Monopole currents are defined in the maximum 
Abelian projection and only one type is shown by 
lines. The upper 3D plot includes the local chiral 
condensate ipiplx), indicated by light grey dots 
whenever a certain threshold is exceeded. One 
clearly sees that both the instantons and antiin- 
stantons are surrounded by a cloud of ipip{ x ) > 
Q. The middle 3D plot exhibits the situation 
for the quark charge density i/j'ip(x) indicated by 
light and dark grey dots depending on the sign 
of the net color charge excess. One observes that 
tp^tp(x) alternates in sign already in one instanton 
implying trivial correlations < tp' ! ip(x)q(y) >= 
(not shown in Fig. 1). The lower 3D plot dis- 
plays the chiral density ipj5tp(x) again indicated 
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Figure 2. Aspects of the fermionic field for a 
fixed time slice of a gauge field with an instanton- 
antiinstanton pair drawn by white and black dots. 
Upper plot: ^t/j(x) > (light grey dots) within 
the instantons. Middle plot: ip^i(;(x) alternates in 
sign (light and dark grey dots) in a single instan- 
ton. Lower plot: tpj^tplx) is positive/negative in 
instanton/antiinstanton. 



by light and dark grey dots. One nicely sees that 
the positive instanton is always surrounded by a 
lump with ip^ipix) > and vice versa. Combin- 
ing the above finding of Fig. 1 showing that the 
correlation functions between fermionic and topo- 
logical quantities are not very sensitive to cooling 
together with the 3D images in Fig. 2, we con- 
clude that instantons go hand in hand with clus- 
ters of 'tjjTtlj(x) ^ , T = 1,74,75, also in the 
uncoolcd QCD vacuum. 

In summary, our calculations of correlation 
functions between topological densities and the 
fermionic observables yield an exponential de- 
crease. Results for the condensate and the mod- 
ulus of the quark charge correlators are almost 
identical, as expected, since the quark condensate 
reflects the absolute value of the quark charge 
density. These correlation functions show little 
cooling dependence. 

The correlations unambiguously demonstrate 
that not only the local chiral condensate but also 
the quark charge and chiral density take non- 
vanishing values predominantly in the regions of 
instantons and monopole loops. Note that for the 
chiral density this behavior is expected due to the 
anomaly. 

Visualization exhibited that the distribution of 
sea-quarks is drastically enhanced around centers 
of nontrivial topology (instantons, monopoles) in 
Euclidian space-time. It must be emphasized that 
this represents the situation on a finite lattice 
with finite quark mass without the extrapolation 
to the thermodynamic and chiral limit. However, 
since all such correlators turned out rather in- 
dependent of the gauge group and choice of the 
action etc., we expect that they are generic. 
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